Exploring correlations in the stochastic Yang— Mills vacuum 
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The correlation lengths of nonperturbative-nonconfining and confining stochastic 
background Yang-Mills fields are obtained by means of a direct analytic path-integral 
evaluation of the Green functions of the so-called one- and two-gluon gluelumps. 
Numerically, these lengths turn out to be in a good agreement with those known from 
the earlier, Hamiltonian, treatment of such Green functions. It is also demonstrated 
that the correlation function of nonperturbative-nonconfining fields decreases with 
the deviation of the path in this correlation function from the straight-line one. 

I. INTRODUCTION 

Stochastic vacuum model [l] (for reviews see 2, 3]) puts in practice the idea that it is the 
stochasticity of nonperturbative background Yang-Mills fields that provides confinement. 
This idea is implemented by assuming the ensemble of fields to be predominantly Gaus- 
sian, with the two-point correlation function containing the part essential for confinement. 
Despite the numerical support by various lattice simulations 3^ 6 ] , this model, appealing 
by its simplicity, requires also theoretical support in the form of analytic studies of the 
said two-point function. Such calculations have been performed in various Abelian models, 
where confinement is provided by the monopole condensation p]], in the instanton- vacuum 
model 8|, and recently also in AdS/QCD [9(. In QCD itself, the first progress has been 
reached in Ref . [kJ , where the correlation length of the two-point function has been related 



to the gluon condensate. Furthermore, Ref. 



Ill ] explored the possibility for nonperturbative- 



nonconfining and confining background fields to have different correlation lengths. Such a 
possibility is specific for QCD, and unlikely to exist in Abelian models with confinement. 
This phenomenon has first been noticed in the lattice QCD simu 



further analytic support within the theory of so-called gluelumps 11 



13|. 



ations |6|J, gaming its 



Gluelumps can be thought as bound states of gluons in the field of a hypothetical infinitely 



heavy adjoint source [ijj]. In Yang-Mills theory, they define the correlation lengths of the 
field-strengths' two-point function in the same way as, in full QCD, physically existing heavy- 



ight 
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mesons define the correlation length of a nonlocal quark condensate (tp(x)Q xx iip(x')) 



15] (here <& xx i is a phase factor along the straight-line path, which is provided by 
the heavy-quark propagator). Unlike the fundamental case, the adjoint case allows for two 
different types of heavy-light objects — those with a single gluon, called one-gluon gluelumps, 
and those with two gluons, called two-gluon gluelumps (cf. Fig. [1] below). The first case 
is similar to the above-mentioned nonlocal quark condensate, whereas the second case is 
conceptually different, as it corresponds to two gluons connected together with the heavy 



source by three fundamental strings. Specifically, it has been shown in Ref. 12] that the 
Green functions of one- and two-gluon gluelumps define respectively the correlation lengths 
of nonperturbative-nonconfining and confining stochastic Yang-Mills fields. Furthermore, 



these Green functions have been explored in Ref. [12] by using respectively one- and two- 
body relativistic Hamiltonians with linear potentials. 

In the next Section, we accomplish the full calculation of quantum-mechanical path inte- 
grals representing those Green functions, which has been started in Ref. [3]. This analysis 
will in particular allow us to find the above-mentioned two vacuum correlation lengths, and 



compare them with those obtained in Ref. 12] within the Hamiltonian approach. Such a 



calculation turns out to be possible by virtue of an effective parametrization of minimal areas 



swept out by the strings in the gluelumps. This parametrization, suggested in Ref. [16], has 
been successfully used there to account for confinement of virtual gluons in the polarization 
operator. In Section III, we summarize the main results of the paper. In Appendix A, we cal- 
culate the Green function of a one-gluon gluelump for paths deviating from the straight-line 
one, and show that the corresponding correlation function of nonperturbative-nonconfining 
background fields decreases with the deformation of the path. 



II. ANALYTIC CALCULATION OF THE VACUUM CORRELATION LENGTHS 

Stochastic vacuum model suggests the following parametrization for the two-point corre- 
lation function of gluonic field strengths 

|^tr (FZ u {x)T a <S> xx ,F b Xp {x')T b <S> x , x ) = (5, x 5 up - S^8 vX )D{x - x') + 
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+ o K (( x ~~ x ')\&pp -{x- x') p 5 vX ) + dl {(x - x') p 8^x -{x- x') x 5 w )] D x (x - x). (1) 
Here, & xx i = V exp \ig f^dz^A^^T^ is a phase factor along some path interconnecting 
the points x' and x. For the rest of the present Section, this path is chosen along the straight 
line. Furthermore, T a 's are the SU(iV c ) generators in the fundamental representation, the 
average (. . .) is taken with respect to the Yang-Mills action ~ J d 4 x(F^ u ) 2 , F^ v = d^A^ — 
d^A^ + g f abc 'A h A c v is the Yang-Mills field-strength tensor, a = 1, . . . , N% — 1. The functions 
D\(x) and D(x) parametrize respectively the nonperturbative-nonconfining and confining 
self- interactions of stochastic background fields. In this Section, we calculate D\(x) and 
D(x) analytically, using their relation to the infra-red Green functions of one- and two- 
gluon gluelumps [121 ]: 

0lW = -4/Q^, DW = ^pi) 5W , (2) 
where C{ = (N% — 1)/ (2iV c ) is the quadratic Casimir operator of the fundamental representa- 

1/2 

tion. Equations ([2]) hold at large distances, |x| > 0{a { ), where cif is the string tension in 



the fundamental representation 



12J. In this infra-red regime, perturbative contributions are 



negligible, and the Wilson-loop averages appearin g in the Green functions G(x) and S(x) 
exhibit the minimal-area law (cf. below). In Ref. 12], the leading infra-red asymptotes of 
these Green functions have been obtained in terms of the lowest eigenvalues of respectively 
one- and two-body Hamiltonians with the linear potential. Below, the full infra-red Green 
functions G(x) and S(x) will be found by means of a direct calculation of the corresponding 
path integrals with the minimal-area ansatze. 

We start with the Green function of a one-gluon gluelump, 

G(x) = J ds J Vz^exp I- J -^d\ - aS min ) , (3) 

where the minimal surface of area S m i n is swept out by the adjoint string of tension a, which 
connects the gluon to the adjointly charged source. While the effects of non-staticity of the 
source are considered in Appendix A, here we assume the source to be static, i.e. evolving 
entirely along the 24-axis. This means that only the ^-coordinate of the point x in Eq. ([3]) 
is nonvanishing, i.e. 

x = (0,L). 

In what follows, we continue with the notations G{x) and S(x), expressing these functions 
in terms of L = \x\. 



Figure 1: A one-gluon gluelump (left) and a two-gluon gluelump (right). For illustration, at various 
values of are also shown a straight-line adjoint string connecting the gluon to the static adjointly 
charged source (for the one-gluon gluelump) and three fundamental strings connecting two gluons 
and the static source (for the two-gluon gluelump). 

We calculate the path integral in Eq. (jHJ) by approximating S m i n through the Cauchy- 



Schwarz inequality [161 ]: 



S min = £dr\z(r)\ < JL^drtf. (4) 

It reduces the path integral to that of a harmonic oscillator with a variable frequency. Indeed, 
using the formula e - ^ = / °° exp (—A — ^) with A > 0, and changing the proper-time 
variable as s — > /j = we have 

where §T>z denotes the integration over trajectories, which start and end up at z = 0. 

az 2 

The kinetic term means, of course, that the auxiliary parameter /x can be viewed as an 
effective gluon mass. Furthermore, the approximate equality "~" in Eq. (jSJ) is understood in 
the sense of approaching the upper limit for 5 m i n in Eq. PJ. We see then that the resulting 
path integral over z±(t) is that of a free particle, 
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whereas the integral over z(r) is that of a harmonic oscillator. It reads 



Vzexp ( — 



drz 2 



a 2 L 

1a" 



drz' 



to 



An sinh ( ^ . } 



3/2 



(6) 



where u = a^j2^iL/\ is the frequency of the oscillator. Changing further the integration 
variable from A to £ = crL 3 / 2 / \ / 2fi\, we can perform the /i-integration: 

a 2 L 3 ' 



G{x) 



a 



L 

32^ 



dfx 



d£ 



ex P I — ?r 



C7 



4vr 2 i V£sinh 3/2 £ 



e where I = aL 2 . 



(7) 



In the infra-red limit Z ^> 1 of interest, this integral can be evaluated analytically as follows: 



% e -i/e + 2 3 / 2 / i^ e -f 4 ~ V + 2 3 / 2 / r 



.Mi 



-/ 



6 A 71 - 



6ctL 



(8) 



Here, the replacement of the lower limit of integration in the second integral by was 
legitimate, since the saddle-point value £ s p = a/2//3, in the limit I ^> 1 at issue, is larger 
than 1 (and therefore the contribution of the integration region < £ < 1 to the whole 
integral is exponentially small). Thus, the leading result in the large- \x\ limit reads 

a 



G(x) 



6cr\x\ 



V3 



Ti ' 



By means of Eq. ([2]), it yields the following function Di(x): 



g 2 a 



2a 

71 



3/2 



6<r|x| 



Dx{x) 

V 71 J \ x \ 

This expression can be compared with the one from Ref. |12| . 

M p -M |x| 



D 1 (x)=g 2 Q 



2 n J "0" e 



2tt 



\x\ 



(9) 



(10) 



(11) 



where the value Mo ~ 1.5 GeV was obtained from the Schrodinger equation with the linear 
potential. To this end, we evaluate the adjoint string tension via the so-called Casimir- 
seeing hypothesis tt This hypothesis, supported both by .att.ee simulations Q aud 
analytic studies [3|, [19[, suggests proportionality of the string tension, in a given represen- 
tation of SU(iV c ), to the quadratic Casimir operator of that representation. For the adjoint 
representation of SU(3), it yields a = |<7f, where the value of the string tension in the 
fundamental representation is <jf = (440 MeV) 2 . Accordingly, the obtained 

mass of the lg gluelump = y/0a ~ 1.6 GeV, (12) 
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turns out to be very close to the quoted value of M . Note that, in principle, a gluon is 
confined up to arbitrarily large distances only at N c ^> 1, whereas at N c ~ 1 it can be 
screened by other gluons. In the large- N c limit, a — > 2<jf, and the exponents in Eqs. ( FlOl) 
and ( 1TT1) coincide numerically extremely well, since in that case y/6a ~ 1.5 GeV. 
We proceed now to the Green function of the two-gluon gluelump: 



S{x) 



dfi f°° dfi 



L 2 



4 Jo /i 2 Jo V 2 



Vz 



x 



o I 1 jo 



drz 2 u - t 

2./0 



-o-f y dr (|z| + |z| + |z - z|)j. 



Trying out the Cauchy-Schwarz inequality, Eq. (j3J), for each of the three distances separately 
would involve integrations over three auxiliary parameters. In order to reduce this number 
to just one, it is more useful to apply the Cauchy-Schwarz inequality in the form 



1 n 
n 



< 



1 n 

. Af, where A { > 0. 

\J i=i 



At n = 3, it yields 



I z I + I z I + I z — z I < Vo • \/z 2 + z 2 + (z — z ) 2 - 



Now, so long as the common square root for the three distances is assembled, we can 
again apply the Cauchy-Schwarz inequality in the form of Eq. (TJl) with only one auxiliary 
integration: 



S(x) 



dfi f°° dfi 



4 Jo V 2 Jo V 2 



Vz„ 



3a 2 L 



Vz„ 



dX 



exp 



-A 



2./,, 



drz 2 u - ± 

2./o 



dTz\- 



4A 



(13) 



dr [z 2 + z 2 + (z - z) 2 ] 

Similar to the one-gluon gluelump, integrations over z±(t) and z±{t) in this formula yield 

-L 



dr{p,z\ + fiz 2 ^ 

We observe now that, if the terms z 2 and z 2 in Eq. (1131) were absent, the path integral 
would be that of two (mutually interacting but otherwise free) particles, which was calculated 





. 2/ 


I X>2 4 / Uz 4 exp 


lo Jo 







(fj, + n)L 


= exp 


2 



in Ref. 



16[. Here, however, we have to deal with two (also mutually interacting) harmonic 



oscillators. Nevertheless, such a path integral § X?z § Dt, can still be calculated. This turns 
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out to be possible upon the diagonalization of the action by virtue of the known fact that 
two positively definite quadratic forms (that are, the kinetic and the potential energies) can 
be diagonalized simultaneously. Passing to the new integration trajectories u(r) and v(r) 
according to the formulae 

z = u + av, z = (3u + v, (14) 
we obtain the diagonalization conditions 

fia + fif3 = 0, a + (3 = (a — l)(/3 — 1). 

The solution to these equations is straightforward: 



/3(/i, u) = —— ■ «(//, u), where a(a, //) = 1 — — + \ — + [ 1 — — ] , (15) 
/' " /' V /' V /'/ 

where we have chosen for concreteness the "+" sign in front of the last square root. Then 
the kinetic- and the potential-energy terms read 

yuz 2 + /2z 2 = (/i + /i/3 2 )u 2 + (/ict 2 + /i)v 2 , 

z 2 + z 2 + (z - z) 2 = 2 [(/3 2 -0+ l)u 2 + (a 2 - a + l)v 2 ] . 

The Jacobian stemming from the change of integration trajectories in Eq. (fl4|) is, of course, 
also /i- and /2-dependent, namely 

VzVz = |1 - a/3| PuDv = ^1 + a 2 -^j VuVv. 

Thus, we have reduced the path integral § T>z § T>z to the product of two well-known path 
integrals for non-interacting harmonic oscillators. Introducing finally the dimensionless vari- 
ables 

- ^ - P A- 



/0f y/(Tf 

we can write down the result in the following form: 



33/2(T f ,5/2 p r dp f, 2»\ r dx 



\2 



cxp 



64vr 9 / 2 J u 3 / 2 J v 3 l 2 V v) Jo A 

x [(/3 2 - /3 + l)(a 2 - a + + z//3 2 )(z7 + z/a 2 )] 3/4 x 



. (v + is)d 
A — 
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Figure 2: Plotted is 



ln(S(x)/af) 



in the range d £ [3, 25]. 



where a = a(i>, v) and = j3(v,i>) are given by Eq. (TT5]) . The remaining ordinary integra- 
tions in this formula have been performed numerically. In Fig. &\ we plot the so-calculated 



HS(x)/af) 



as a function of d in the range d e [3, 25], and observe an asymptotic approach 



of this quantity to 6.0 at large d. For er f = (440 MeV) 2 , the corresponding 



mass of the 2g gluelump = 6y/af ~ 2.6 GeV 



(16) 



agrees remarkably well with its value of 2.56 GeV found in Ref. 12J within the Hamiltonian 
approach. 



III. CONCLUDING REMARKS 



In this paper, we have analytically confirmed the statement that nonperturbative- 
nonconfming and confining self-interactions of stochastic background Yang-Mills fields can 
have different correlation lengths. Namely, we have obtained these lengths from the Green 
functions of one- and two-gluon gluelumps, by means of an analytic treatment of the 
quantum- mechanical path integrals for one and two gluons (inter) connected by strings with 
the static adjoint source. The resulting inverse correlation lengths are given by Eqs. (1X2 ) 
and f|T6|) . They turn out to be in a remarkably good agreement with those found in Ref. [12 1 
by means of a different, Hamiltonian, approach. We have also demonstrated that, in ac- 
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cordance with the lattice results [5], the two-point correlation function of nonperturbative- 
nonconfining background fields decreases with the deformation of the path interconnecting 
these two points. On the technical side, for the calculation of the Green function of a two- 
gluon gluelump, a novel method has been developed of an analytic evaluation of the path 
integral with the common minimal surface formed by two gluons together with the adjointly 
charged source. 

In the forthcoming studies, we plan to incorporate dynamical quarks in this theoretical 
framework. Their appearance can lead to the breaking of not only fundamental strings, but 
also of the adjoint ones — through the creation of q-g-q hybrids. An effective reduction, 
due to the string-breaking, of string world sheets in the spatial directions can lead to their 
enhancement in the temporal direction, that is, to the increase of the vacuum correlation 
lengths. In this way, vacuum correlation lengths in full QCD can exceed those obtained here 
in pure Yang-Mills theory by a factor of 3, as suggested by the lattice simulations ^| and 
the phenomenology of hadronic collisions 20] . 
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Appendix A: Path-dependence of the function D\{x) 

Lattice simulations [5] indicate that, if the path interconnecting the points x and x' 
in Eq. ([T]) deviates from the straight-line one, the correlation function decreases. In this 
Appendix, we test this indication analytically, at an example of the function Di(x), by 
evaluating it for paths of various shapes. The approach we use for this study is the same as 
the one used in Section II. 

We parametrize the path interconnecting x and x' by some vector-function w(r), whose 
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concrete form will be specified below. Accordingly, the minimal-area ansatz fll]) gets modi- 
fied, and takes the form 



S min = £ dr|z(r)-w(r)| < J L £ dr(z — w) 2 . 

The emerging path integral §T>z turns out to be similar to that of Eq. (jHJ), and reads (cf. 
Ref. |2l|) 

3/2 



Dzexp I / drz 



At r , . ., n -7. r/ - 



4A 



J 2 ^ f L , \ 

drz H / ctrzw 



CO 



4tt sinh ( ^ 



x exp 



-; — r- / du dv sinh 

8/i 2 sinh fgj -A) Jo 



.2/1 



sinh ( ^— f ) w(Ltt)w(Lu) ^ . 



The corresponding path-dependent Green function G w (x) generalizes Green function G(x 
and can be written as [cf. Eq. ([7])] 

a 2 L 3 



G w (x) = a 



The function 



l r°° dp r d A / ^M4.f^i 



(Al) 



— - — I du f dv sinh[£(l — u)] sinh(^v)w(Lu)w(Lv) 
smh^Jo Jo 



(A2) 



introduced here, vanishes in the limiting case w = 0. Performing again the /i-integration 
exactly, and using Eq. with G{x) replaced by G w (x), we obtain the path-dependent 
correlation function 



D l '(x)=a — 



7T" 



(^sinhO 372 



e « 



(A3) 



where I is defined in Eq. ((7j). To calculate this integral, we parametrize the path w(r) in 
the form which provides smooth approximations to step-like paths used in Ref. jfj. Namely, 
we consider two types of smooth paths: 



w«(r) = ^(sin^,0,0), w^(r) 



L ( 2vrr 

sin , 0, 



2 V L 

Each of these paths lies in the (l,4)-plane, deviating from the £4- axis to the maximum 
distance L/2, which is still compatible with L. In general, for larger deviations, one can 
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/*1 r-u 

' du dv sinh[£(l — u)} sinh(£?;) sin(7rw) sin(7rt> ] 
o Jo 



expect cancellations of contributions stemming from the mutually backtracking pieces of the 
path. 

Consider first the path The corresponding function ( 1A2I) reads 

MW 4 \2 sinh£ 

The ^-integration in Eq. (IA3[) can be performed in the same way as in Eq. (JSJ), by splitting the 
integration region into the intervals [0, 1] and (1, oo). In the interval [0, 1], /(£) = ^- + 0(£ 4 ), 
and it can be disregarded altogether compared to 1. The contribution stemming from 
this interval again appears exponentially suppressed compared to the leading contribution 
stemming from the interval (1, oo). For £ G (1, oo), the exponential part of the ^-dependence 
of the integrand, which determines the position of the saddle-point, reads [cf. Eq. 



e 2 £ \f 1+ ftt) _ rp Q gg ure ou ^ foe extent to which the function /(£) affects the saddle point, 
one needs to consider the corresponding saddle-point equation, 

---^v/l + /(0 + - pQS = = 0, (A4) 

2 £ 2 2£ V /TT7(0 



where 



4£e ^ / efot / dv sinh[£(l — u)} sinh(£t> ) sin(7ra) sin(7r?;) 

2 Jo Jo 



at £ > 1. (A5) 



One can then prove numerically that the absolute value of the second term on the left-hand 
side of Eq. ( 1A4I) exceeds the third term by at least one order of magnitude. Therefore, 
the third term can be disregarded compared to the second one, for the reason of smallness 
of /'(£)• Accordingly, approximating /(£) by some constant /, we obtain the saddle-point 



value £ s . p . = ^J2ly/T~+ //3, which is much larger than 1 in the limit \fl ^> 1 of interest. 
And indeed, for £ ^> 1, one can see a very weak variation of the function /(£), which is 
illustrated by Fig. El where this function is plotted in the range £ G [5,30]. We therefore 
approximate /(£) by its limiting value at £ ^> 1, i.e. set / ~ 1.23. 

Then, given that £ s p . lies deeply inside the region (1, oo), we can replace the lower limit 
of integration by [cf. Eq. (jSJ)], and obtain 

Comparing this expression with Eq. (flOj) . we observe a decrease of the amplitude of the 



correlation function by a factor of + f ~ 1.22, in a qualitative agreement with indications 
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1.25 




Figure 3: Function f|A5j) for £ » 1. 

of Ref. [5] on the dominance of the straight-line path. Furthermore, we obtain also an 
increase of the mass of the one-gluon gluelump, Eq. ffl2|) . by the same factor of 1.22. Such 
an increase yields an even stronger suppression of contributions produced by curved paths. 

In general, the present approach yields, for a given path, a decrease of the vacuum 
correlation length and the amplitude of the correlation function by the same amount as 
compared to these quantities for the straight-line path. With the deformation of the path, 
the suppression factor increases. This statement can be illustrated by considering the path 
w( 2 )(r), which is deformed stronger than w^(r). The corresponding function 



/(0 - \ 



du I dv sinh[£(l — u)\ sinh(£t> ) sin(27ru) sin(27rt> ) 



at £ > 1 



also exhibits a rapid levelling-off, similar to its counterpart given by Eq. ( 1A5I) [cf. Fig. |3]. 
The difference is that the limiting value of the function /(£) at £ ^> 1 appears in this case 
larger, namely / ~ 4.86, instead of 1.23. Accordingly, the suppression factor in Eq. f ]A6j) 
becomes \/l + f ~ 1.56, instead of 1.22. Such an increase of the factor yfl + f is quite 
fast, in spite of its slow, fourth-root, functional dependence. Thus, our analysis suggests a 
substantial suppression of contributions to the two-point correlation function, which stem 
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from strongly deformed paths. 
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